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ABSTRACT 

We have studied the effect of adopting different values of the total baryonic mass surface density in 
the local disk at the present time, £(i? , teal), m a model for the chemical evolution of the Galaxy. We 
have compared our model results with the G-dwarf metallicity distribution, the amounts of gas, stars, 
stellar remnants, infall rate and SN rate in the solar vicinity, and with the radial abundance gradients 
and gas distribution in the disk. This comparison strongly suggests that the value of £(-Rq, few) which 
best fits the observational properties should lie in the range 50-75 M pc -2 , and that values of the 
total disk mass surface density outside this range should be ruled out. 

Subject headings: Galaxy: evolution — Galaxy: stellar content — Galaxy: structure — dark matter 



1. INTRODUCTION 



Despite a number of estimates, obtained by means of 
several different methods, the value of the local column 
density of the Galactic disk in stars and gas, Sb(i?© ~ 



8 kpc) = T, sta 



-'gas j 



has been until recently very un- 



certain. Direct HST star counts imply a stellar column 
density of ~ 27 M pc~ 2 (Flynn, Gould, & Bahcall 1999), 
while the (H I + H2) column density amounts to 7-13 
M Q pc- 2 (Dickey 1993; Flynn et al. 1999). Therefore, 
the surface density of the directly seen baryonic matter 
amounts to 35-40 M pc -2 . However, the actual column 
density of the "luminous" disk could be higher because the 
above detection does not account for dead stars and brown 
dwarfs. 

The total disk surface density E z (i? ), including also 
any eventual disk dark matter component, can be deter- 
mined from the motion of stars in the solar neighbor- 
hood perpendicular to the Galactic plane (e.g. Binney 
& Tremaine 1987; Gilmore, Wyse, & Kuijken 1989). From 
the kinematics of a sample of K-giants Bahcall, Flynn, & 



Gould (1992) derived E Z (R(. 



84 



+30 
-25 



M Q pc" 2 (at 1 er 



level) . Notice that such a measure is a complex one and it 
is strongly dependent on the underlying assumptions. By 
employing different samples and/or analysis, lower values 
for E z (-R ) have been claimed: £ z (i? ) = 46 ± 9 M 
pc~ 2 (Kuijken & Gilmore 1989) — revised to 48 ± 9 M 
pc" 2 (Kuijken & Gilmore 1991), E z (i? ) = 52 ± 13 M 
pc -2 (Flynn & Fuchs 1994). The local column density 
of stars and gas could be smaller, because these dynami- 
cal estimates may include some non-baryonic matter. No- 
tice however that it is quite unlikely that there is non- 
baryonic matter in the disk, because non-baryonic dark 
matter is probably non dissipational. Any non-baryonic 
matter present in the disk is likely to be part of the dark 
halo, and this component is always subtracted from the 



dynamical estimates of the local column density. 

Fitting the Galaxy rotation curve (RC) with a mass 
model provides a further way of estimating the local col- 
umn density of disk matter Ey. In fact, such a fit indi- 
cates what is the fraction /3 2 of the circular velocity due 
to the disk component (at the Sun's position). Then, 
V 2 is compared with the circular velocity of an expo- 
nential thin disk, Vj sfc = G^Ey Rd f(R Q /R D ), with / 
a known function of the Galactocentric distance of the 
Sun expressed in terms of disk length scales (see Freeman 
1970), to yield Ey, providing that i? and Rrj are known. 
The results are far from unique, and reflect the presence 
of uncertainties in the mass modelling, in the actual shape 
of the Galactic RC, and in the basic structural parame- 
ters. The values found range between Ey ~ 50 M pc -2 
and - 100 M pc~ 2 (Oiling & Merrifield 1998; Dehnen 
& Binney 1998), and moreover include any dark matter 
distributed like the stellar disk. In external galaxies of lu- 
minosity similar to the Galaxy (Lb ~ 10 11 Lq), values of 
Ey as high as ~ 90 M pc -2 at R = 8 kpc cannot be 
ruled out (Persic, Salucci, & Stel 1996). 

Finally, let us note that if the dark matter halo in the 
Galaxy follows the standard cold dark matter (CDM) uni- 
versal density profile (Navarro, Frcnk, & White 1997), and 
therefore at the Sun position pscdm oc R~ 2 , then the lo- 
cal disk column density is required to be towards the low 
end of the above-cited values (e.g. Courteau & Rix 1999). 
On the other hand, low-fi CDM density profiles are not in- 
consistent with a baryon dominated region inside i? , and 
then with much higher values for the local disk column 
density (Navarro 1998). 

Summarizing: all the above observations, often entan- 
gled with heavy theoretical assumptions, poorly constrain 
the value of the column density of the baryonic matter in 



^^Dipartimento di Astronomia, Universita di Trieste, via G.B. Tiepolo 11, 1-34131 Trieste, Italy; matteucci@ts.astro.it 
2 It should be noted that the value for /3 depends on the assumed profiles for the disk and the halo. 
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the local disk, that may lie between 35 and 100 Mq pc . 

However, very recently the results of the European 
Space Agency's Hipparcos mission have allowed a drasti- 
cal reduction of this permitted range, by making available 
direct distances for the tracer stars. In particular, Creze 
et al. (1998) and Holmberg & Flynn (1998) have analyzed 
the A and F stars in the Hipparcos data set (some 10,000 
stars) and derived greatly improved estimates of the total 
gravitating mass: the gravitating disk mass seems to be 
now firmly established at 50 to 60 Mq pc~ 2 . 

In this paper we investigate how the gravitating mass of 
the disk influences the Galactic chemical evolution. Suc- 
cessful models for the chemical evolution of the Galaxy, 
in fact, require a star formation law depending on the to- 
tal mass surface density (Tosi 1988; Matteucci & Francois 
1989; Chiappini, Matteucci, & Gratton 1997), and the star 
formation rate (SFR) is a fundamental parameter in the 
evolution of galaxies. 

2. THE CHEMICAL EVOLUTION MODEL 

We adopt the two-infall chemical evolution model of 
Chiappini et al. (1997), to which we refer for an exhaustive 
explanation of the main assumptions and basic equations. 
Here we just review the features of the model most directly 
related to the aspects discussed in the present study. 

Briefly, the overall evolutionary scenario is the follow- 
ing: it is assumed that the Galaxy formed out of two main 
infall episodes; during the first episode, the primordial gas 
collapsed very quickly and formed the spheroidal compo- 
nents (halo and bulge); during the second episode, the 
thin disk formed almost independently of the previous in- 
fall episode, mainly out of material of primordial chemical 
composition. The disk formation process is assumed to 
be slower with increasing Galactocentric distance (Larson 
1976; Matteucci & Francois 1989); this "inside-out" forma- 
tion of the Galactic disk is required to reproduce the abun- 
dance gradients and the gas distribution along the disk. A 
hiatus in the star formation between thick disk and thin 
disk phases, as suggested by recent observations (Grat- 
ton et al. 1996; Bernkopf & Fuhrmann 1998; Fuhrmann 
1998; Carraro 2000), is naturally produced by the model 
under the assumption that the star formation stops when 
the gas surface density drops below a threshold of 7 Mq 
pc -2 (Kennicutt 1989). Moreover, the presence of such a 
threshold predicts a roughly constant gas density profile 
at the outer radii (between 8 and 14 kpc) in agreement 
with observations (Dame 1993). We will focus on the thin 
disk formation process. 

The most relevant quantities in the present discussion 
are the infall rate and the SFR, since they involve the total 
mass surface density profile. The rate at which the thin 
disk is formed out of external matter (although it has also 
some initial contribution from the halo gas) is: 

glM = B(fi)e-("-)/ T °, (1) 
at 

where S/(i?, t) is the gas surface density of the infalling 
material, which has a primordial chemical composition; 
tmax is the time of maximum gas accretion onto the disk, 
coincident with the end of the halo-thick disk phase which 
is set equal to 1 Gyr; m is the timescale for mass accre- 
tion onto the thin disk component. We assume that td 



increases with increasing the Galactic radius: 

t d (R) = 1.033 x R- 1.267. (2) 

This relation derives from the fact that in order to fit the 
G-dwarf metallicity distribution in the solar neighborhood 
it is necessary to assume tjj = 7 Gyr, whereas to repro- 
duce the metallicity distribution of the stars in the Galac- 
tic bulge (R = 2 kpc) a timescale td — 0.8 Gyr is required 
(Matteucci, Romano, & Molaro 1999). It should be em- 
phasized that this choice also guarantees a good fit to the 
gas density profile and to the Galactic abundance gradi- 
ents. The quantity B(R) is derived from the condition of 
reproducing the current total mass surface density distri- 
bution along the disk: 

B(R) = , ^ tG ^ ( 3 ) 

T D (R)(l - e-(*Oa|-*m..)/T D (fl)j 

teal = 15 Gyr is the age of the Galaxy starting with the 
formation of the halo. From equation (3) it is clear the 
strong dependence of the infall rate from the total mass 
density profile. 

Observations in our own and in other disk galaxies al- 
low us to relate the SFR to intrinsic parameters of galax- 
ies, in particular to the average surface gas density: SFR 
cx ~Eg as (Schmidt 1959), where the exponent is reason- 
ably well known (e.g. k = 1.4 ± 0.15 — Kennicutt 1998). 
Moreover, the correlation between metallicity and surface 
brightness noted for late type spirals (McCall 1982; Ed- 
munds & Pagel 1984; Dopita & Ryder 1994) is consistent 
with theories of self-regulated star formation, in which the 
energy produced by young stars and SNe feeds back into 
the ISM and inhibits further star formation by producing 
an expansion of the region surrounding the new stars or the 
SN remnants, thus relating the SFR to the gravitational 
potential and therefore to the total mass surface density. 
In addition, according to the most popular prescriptions, 
one of the key conditions for star formation to occur is 

t cooling ^ffi (^) 

where t coo u ng is the cooling timescale and tff is the free- 
fall timescale (Katz 1992; Navarro & White 1993; Carraro, 
Lia, & Chiosi 1998); tff oc p^ 1 ^ 2 , where p is the total mass 
density (Buonomo et al. 1999). Again, in some way the 
SFR is related to the potential well and therefore to the 
total mass surface density. The SFR adopted here is the 
same as in Chiappini et al. (1997) and it was chosen in 
order to give the best agreement with the observed con- 
straints: 

t). (5) 

v(t) is the efficiency of the star formation process, which is 
set to 1 Gyr -1 , except when the gas surface density drops 
below 7 M Q pc~ 2 . In fact in this case v{t) = 0, because 
below this critical density the gas is gravitationally stable 
against density condensations into stars. t) is the to- 
tal disk mass surface density at a given radius R and a 
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Table 1 

Comparison between model predictions and observations for some relevant current quantities. 



Quantity 


35 


54 


80 


100 


Observations 


Reference 


Type I SNe (century -1 ) 


0.27 


0.37 


0.59 


0.72 


0.17 


-0.7 


1 


Type II SNe (century -1 ) 


0.73 


0.95 


1.67 


2.06 


0.55 


-2.2 


1 


Nova outbursts (yr -1 ) 
SFR{R Q ,t Ga i) (Ms pc- 2 Gyr- 1 ) 


20 


27 


35 


41 


20- 


-30 


2 


0-3 


2.57 


2.89 


3.58 


2- 


10 


3 


T, gas (R Q ,t G al) (Mg pc^ 2 ) 


7.0 


7.0 


8.4 


10.3 


7- 


13 


4 


S siars (R Q ,t Ga i) (M s pc" 2 ) 


25.5 


38.4 


54.3 


65.6 


35: 


1=5 


5 


'Eremnants (R®,t Ga l) (Mg pC ) 

t infa ii{RQ,t Gal ) (M Q pc~ 2 Gyr" 1 ) 


2.5 


8.6 


17.3 


24.0 


2-4 


6 


0.57 


0.89 


1.32 


1.66 


0.3- 


-1.5 


7 



References. — (1) Chiappini et al. 1997 and references therein; (2) Shafter 1997; (3) Glisten & Mezger 1982; (4) Dickey 1993; Flynn et al. 
1999; (5) Gilmore et al. 1989; (6) Mera, Chabrier, & Schacffer 1998 quote T,wDs+NSs = 2-4 Mq pc -2 ; here S remnan t s = 'EwDs+NSs+BHs', 
(7) Portinari, Chiosi, & Bressan 1998 and references therein 



given time t, £(i? ,i) is the total disk mass surface den- 
sity at the Solar position and T, gas (R, t) is the gas surface 
density. Note that the gas surface density exponent, k, 
equal to 1.5, was obtained from the best model of Chi- 
appini et al. (1997) in order to ensure a good fit to the 
observational constraints at the solar vicinity. This value 
is in good agreement with the recent observational results 
of Kennicutt (1998) and with N-body simulation results 
by Gerritscn & Icke (1997). 

As far as our present model is concerned, we stress that 
the SFR along the Galactic disk, and hence the evolution- 
ary history of the Galactic disk itself, is fixed by the specific 
value of the local column density of disk matter that we 
adopt. In fact, the total disk mass surface density profile, 
an exponential with scale length Rd, can be expressed as: 

X(R,t Ga i) = X(Re,t G ai)-e-( R - R °y RD , (6) 

which relates the column density at a radius R to the col- 
umn density at the Solar position. For the purpose of this 
paper, we let the local mass surface density value to vary 
from 35 to 100 M© pc~ 2 and we set Rq — 8 kpc. 

Finally, the stellar nucleosynthesis prescriptions we 
adopt are from: 

- van den Hoek & Groenewegen (1997) for low- 
intermediate mass stars (0.8-8 M ) (their case with 
variable ri AGB ); 

- Charbonnel & do Nascimento (1998) for 3 He produc- 
tion/destruction in low-mass stars (M < 2 M ); 

- Woosley & Weaver (1995) for Type II SNe (M > 10 

M ); 

- Thielemann, Nomoto, & Hashimoto (1993) for Type 
la SNe (exploding white dwarfs in binary systems as 
described in Matteucci & Greggio 1986); 

- Jose & Hernanz (1998) for classical novae. 

Our prescriptions differ from those of Chiappini et 
al. (1997) in: i) the range of low- intermediate masses (0.8- 
8 M Q ), where they adopted the Renzini & Voli (1981) 

•^Ionization corrections for low-metallicity H II regions suggest that 
Grucnwald, & Stcigman 1999). 



yields; ii) the fact that we are including the explosive nu- 
cleosynthesis from nova outbursts (Romano et al. 1999; 
Romano & Matteucci 2000); iii) using extra-mixing pre- 
scriptions for the 3 He nucleosynthesis (Chiappini & Mat- 
teucci 2000) and iv) the Sun position, now taken as 8 kpc 
instead of 10 kpc. 

As far as the composition of the infalling material is con- 
cerned, we set Yp = 0.245 rather than 0.23 (where Yp is 
the primordial mass fraction in form of 4 He) , according to 
a recent estimate of the pristine 4 He abundance obtained 
as a weighted mean of the He mass fractions from the two 
most metal-deficient blue compact galaxies, corrected for 
the stellar contribution (Y P = 0.245 ± 0.002 — Izotov et 
al. 1999 3 ). The primordial abundances by mass of deu- 
terium and helium-3 are Dp = 4.5 x 10~ 5 and 3 Hep = 4.0 
x 10~ 5 , respectively. The adopted IMF is that of Scalo 
(1986). 

3. MODEL PREDICTIONS 

3.1. The Solar Circle 

In Tables 1 and 2 the model predictions together with 
the corresponding observed quantities are shown. Differ- 
ent models have identical parameters apart from the value 
ofS(i? ). _ 

From an inspection of Table 1, we see that there are not 
big differences among the quantities predicted by the dif- 
ferent models, with the exception of the surface densities 
of live stars and stellar remnants, which seem to exclude 
the models with the highest total disk mass surface densi- 
ties (80-100 M Q pc~ 2 ). The lowest density value (35 M 
pc~ 2 ) predicts a stellar density which is only in marginal 
agreement with the observed one. 

The predictions on the abundances of the interstellar 
medium at the time of the Sun's formation are indepen- 
dent of the chosen total disk mass surface density but de- 
pend on all other model assumptions, in particular on the 
assumed stellar yields. Therefore, since we adopted a set 
of new yields relative to Chiappini et al. (1997), a compar- 
ison between our results for S = 54 M pc~ 2 and those 

this estimate should be reduced to Yp = 0.241 ± 0.002 (Viegas, 
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Table 2 

Theoretical solar abundances (mass fraction) are compared to the observed ones (Anders & Grevesse 1989). 



Xi 


This paper 


Phinrmini pt al 1 QQ7 


Observed 


H 


0.710 




0.731 


0.706 




13 


3 931 x 10" 

O.^tJ-L A -L\J 


5 


4 630 x 10 -5 


4 80 x 10" 


5 


3 He 

Ale 


3 977 x 10" 


5 


10 01 v 1 -5 


9 93 x 10" 


5 


4 He 


0.273 




0.255 


0.275 




7 Li 


8 740 x 1 0" 


9 




9 63 x 10" 


9 


12 C 


3 Q20 x 10" 


3 


1 827 x 10 -3 


3 03 x 10" 

fJ'KJtJ /n -L\J 


3 


13 c 


5.215 x 10" 


5 


4 758 x 10 -5 


3.65 x 10" 


5 


14 N 


1.670 x 10" 


3 


1.386 x 10~ 3 


1.11 x 10" 


3 


16 Q 


7.378 x 10" 


3 


7.278 x 10~ 3 


9.59 x 10" 


3 


20 Ne 


9.729 x 10" 


4 


9.42 x 10~ 4 


1.62 x 10" 


3 


24 Mg 


2.563 x 10" 


4 


2.48 x 10" 4 


5.15 x 10" 


4 


Si 


7.348 x 10" 


4 


7.03 x 10 -4 


7.11 x 10" 


4 


s 


3.220 x 10" 


4 


3.071 x 10" 4 


4.18 x 10" 


4 


Ca 


4.145 x 10" 


5 


3.95 x 10" 5 


6.20 x 10" 


5 


Fc 


1.467 x 10" 


3 


1.37 x 10~ 3 


1.27 x 10" 


3 


Cu 


8.920 x 10" 


7 


8.18 x 10" 7 


8.40 x 10" 


7 


Zn 


2.635 x 10" 


(i 


2.44 x 10~ 6 


2.09 x 10" 


6 


z 


1.700 x 10" 


2 


1.433 x 10" 2 


1.886 x 10 


-2 



of Chiappini et al. (1997) is required. Table 2 shows that 
the largest differences in the predicted solar abundances 
arise for D, 3 He, 12 C, and 13 C. The improvement in the 
predicted solar abundance for 12 C and 13 C is due to the 
fact that we adopt the new nucleosynthetic yields from van 
den Hoek & Groenewegen (1997) for low and intermedi- 
ate mass stars, whereas the improvement in the predicted 
solar abundances for D and 3 He is due to the fact that we 
take into account that 93% of evolved stars undergo the 
extra-mixing on the RGB and thus destroy, at least partly, 
their 3 He (Charbonnel & do Nascimento 1998). 

One of the most severe observational constraints to 
chemical evolution models is the relative number of stars 
formed at a given time, i.e. at a given metallicity. In Fig.l 
we compare the predicted and observed G-dwarf metal- 
licity distributions; in particular, we show the results of 
models in which £(i?o, tc a i) is taken to be 20, 54 and 100 
M G pc" 2 , respectively. Models with £ = 35 and 80 M 
pc -2 have also been computed. They produce theoretical 
G-dwarf distributions very similar to those obtained with 
S = 54 or 100 Mq pc" 2 , so we do not show them here. 
From an inspection of Fig.l we see that the case E = 20 
Mq pc -2 is in trouble in reproducing the high metallicity 
tail of the distribution while overestimating the number of 
low metallicity stars, whereas all the cases with E lying 
between 35 and 100 Mq pc" 2 lead to a better agreement 
with the observations. We performed a x 2 test and ob- 
tained the following results: \ 2 = 0-H f° r t ne case E = 
20 M G pc" 2 , x 2 = 0.04 for the cases E = 54 and 100 Mq 
pc -2 . We conclude that the observed differential metallic- 
ity distribution is not a useful tool to distinguish among 
realistic values for the total baryonic mass density in the 
local disk. This is not surprising since the G-dwarf metal- 
licity distribution depends mostly on the assumed law for 
the disk formation (i.e. on the assumed infall timcscale). 

Fig. 2 illustrates the different temporal behaviour of 
the SFR under different choices of E. In the low den- 
sity case a stochastic SFR (between and 3 Mq pc -2 



Gyr -1 at the present time) is expected during most of 
the evolution, due to the fact that the threshold in the 
gas density, below which the star formation stops, is eas- 
ily attained in this case. On the other hand, in the 
case of a very high E the star formation process is not 
able to consume enough gas to go below the thresh- 
old. This explains why the highest E values produce 
the high stellar and remnant densities listed in Table 1. 
In a recent paper Rocha-Pinto et al. (2000) suggest that 




-1 -0.5 

[i'c/H] 

Fig. 1. — Theoretical G-dwarf distributions obtained by choosing 
different values of T,(R,Q,tQ a [) as input parameters in the chemical 
evolution model compared to the observed one (Rocha-Pinto & Ma- 
ciel 1996). The [Fe/H] ratios are normalised to the theoretical solar 
one for each model. 



5 




5 10 15 



Time (Gyr) 

Fig. 2. — Temporal evolution of the star formation rate in the 
solar vicinity as predicted by models adopting T,(RQ,tQ a i) = 35, 54, 
and 100 M pc" 2 . 



the SFR in the disk should have proceeded in three dis- 
tinct bursts at 0-1, 2-5, and 7-9 Gyr ago. The ear- 
liest episode of enhanced star formation should reflect 
the beginning of the solar vicinity formation (and prob- 
ably corresponds to our maximum infall time). The other 
two episodes could be related to cyclical mechanisms (for 
instance, tidal interactions between the Galaxy and the 
Magellanic Clouds, or even induced star formation by spi- 
ral arms). As already explained, the theoretical G-dwarf 
metallicity distribution is mainly dependent on the as- 
sumed timescale for the formation of the solar vicinity, 
hence it would not be too much affected by recent bursts 
of star formation. 

The elemental abundance ratios as a function of metal- 
licity are nearly independent of the adopted S value, de- 
pending mostly on the nucleosynthetic yields and the ini- 
tial mass function. Therefore we do not discuss them here. 



3.2. Radial Profiles 

Here we discuss the predicted radial properties, in or- 
der to investigate the model behaviour along the overall 
Galactic disk under different prescriptions on the value of 
EOR ). 

We assume that the radial surface density distribu- 
tion of the disk is exponential (see Sect. 2, equation (6)), 
and run models assuming different values for the expo- 
nential scale length, Rjj. Fig. 3 shows the radial dis- 
tribution of the present day surface gas density as pre- 
dicted by models adopting Rd = 2.5, 3.5, and 4.5 kpc 
(dashed, solid, and dotted lines, respectively), for differ- 
ent S values. The model predictions are compared with 
the observed gas distribution derived by Dame (1993). 

4 A scale length value of 3.5 kpc is the preferred one, being able to 
one (2.5-3.0 kpc — Sackctt 1997; Freudenreich 1998). 




5 10 15 



R (Kpc) 

Fig. 3. — Radial distribution of the present day surface gas den- 
sity as predicted by models with Kr> = 2.5 (upper panel), ~Rr> = 
3.5 (middle panel), and Kr> = 4.5 (lower panel). From bottom to 
top, for each of these panels: lines referring to T,(Rq, ta a l) = 35, 54, 
80, and 100 Mq pc (dotted, continuous, long-dashed, and short- 
dashed lines, respectively). The observed (H I + H2) distribution 
is taken from Dame (1993). The surface density distribution of the 
total gas Tigas is obtained from the sum of the H I and H2 distri- 
butions, Eijj + Ehoi accounting for the helium and heavy elements 
fractions (thick line). 



All our models predict a flattening of the distribution at 
radii larger than 10 kpc, due to the presence of the thresh- 
old star formation density. Models with £ = 35 Mq pc -2 
give a flat gas distribution already at the smallest radii; on 
the contrary, models with S = 100 M Q pc -2 produce very 
steep distributions, hardly compatible with observations. 
All models predict that the gas content increases toward 
the Galactic center, whereas the observations clearly show 
a decreasing trend below 5 kpc. This problem could be 
solved by including the Galactic bar (e.g. Portinari & 
Chiosi 2000), not present in the model. Taking into ac- 
count these uncertainties we can conclude that the best 
gas distribution is obtained for local values of the total 
mass surface density in the range 50-75 M Q pc~ 2 . 

We also notice how increasing the disk scale length from 
2.5 to 4.5 leads to a flattening in the predicted distribution. 
This behaviour is also observed in the case of the predicted 
radial density distributions of stars 4 , stellar remnants, and 
in the case of the predicted abundance gradients. In par- 
ticular, higher total mass densities and larger scale lengths 
tend to flatten the gradients in the external regions of the 
disk. The actual situation is still controversial, with ob- 
servations partly supporting this flattening (Vflchcz & Es- 
teban 1996; Maciel & Quireza 1999) and partly not (Af- 
flcrbach, Churchwell, & Werner 1997; Smartt & Rolleston 
1997; Rudolph et al. 1997). Therefore we can not draw 
any firm conclusion on this point. 

:uce a final stellar density scale length in agreement with the observed 
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4. CONCLUSIONS 

We used a detailed model of galactic chemical evolution 
applied to our own galaxy in order to put constraints on 
the total amount of disk baryonic matter. Observational 
estimates of this quantity in the local disk span a wide 
range of values (from 35 to 100 M Q pc~ 2 ), depending on 
the underlying theoretical assumptions and on the adopted 
methods. We show that by means of chemical evolution 
models it is possible to substantially restrict the observed 
range of the total mass surface density. 

We show that the value of £(i?©) should be restricted to 
the range 50-75 M Q pc~ 2 , in order to guarantee the best 
fit to all the observational constraints available for the so- 
lar neighborhood and the overall Galactic disk. This is well 
consistent with the mass range which Creze et al. (1998) 



and Holmberg & Flynn (1998) advocate from their anal- 
ysis of A and F stars with parallaxes and proper motions 
from the Hipparcos satellite. 
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